Dynamical effects of a cosmological constant 
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Abstract 

The observational evidence for the existence of a non-zero cosmo- 
logical constant is getting stronger. It is therefore timely to address 
the question of its eventual effect on the dynamics of galaxies, clus- 
ters and larger structures in the Universe. We find, contrary to a 
recent claim, that the influence of the cosmological constant has to be 
negligible for, e.g., the rotation curves of galaxies. On larger scales, 
starting with large galaxy clusters, there are potentially measurable ef- 
fects from the repulsive addition to the Newtonian gravitational force 
caused by the cosmological constant. 



During the past few years, remarkable progress has been made in cosmol- 
ogy, both observational and theoretical. One of the outcomes of these rapid 
developments is the increased confidence that most of the energy density of 
the observable universe is of an unusual form, i.e., not made up of the ordi- 
nary matter (baryons and electrons) that we see around us in our everyday 
world. 

There are convincing arguments for the existence of a large amount of 
non- luminous, i.e., dark, matter. The matter content of the universe is at 
least a factor of 5 higher than the maximum amount of baryonic matter 
implied by big bang nucleosynthesis. This dark matter is thus highly likely 
to be "exotic", i.e, non-baryonic. 

There are also indications, although still not entirely conclusive, of the 
existence of vacuum energy, corresponding to the famous "cosmological con- 
stant" that Einstein introduced but later rejected (although without very 
good reasons) in his theory of general relativity. This possibility has recently 
been given increased attention due to results from Type la supernova surveys 
]!], 0. (For a recent review of the observational status of dark matter and 
dark energy, see ||.) 

It may be interesting to investigate possible consequences of the cosmo- 
logical constant besides its influence on the geometry of the universe, and the 
redshift- dependence of the luminosity-distance relation for standard candles 
This is the subject of the present paper. We find disagreement with a 
recent paper [H where the cosmological constant was claimed to influence 
the rotation curves of galaxies strongly. However, some small effects on the 
dynamcis of galaxy clusters do not seem excluded. 

Let us first set our conventions. Einstein's equations read 



The energy density in the form of a cosmological constant A can be conve- 
niently written in units of the density scaled to the critical density, 
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with Gn = l/?7ipi (the numerical value of the Planck mass is mp\ = 1.2 • 10 19 
GeV) and the present value of the critical density 

Thus, 

= Jk- (5) 
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Writing H = lOO/i kms 1 Mpc 1 (with h ~ 0.6 ± 0.1 from observations), 
the present numerical value of the critical density is 



p crit ~ 8 • 10- 47 h 2 GeV 4 . (6) 



This was derived using particle physics units (c = % = 1). Expressed in 
cgs units, the presently observationally favoured value || fi A ~ 0.7 ± 0.2 
then translates to 

A obs ~ 10~ 56 cm" 2 . (7) 

In a recent preprint ||, an attempt was made to explain the flat rotation 
curves of galaxies as being due to the effect of a cosmological constant instead 
of the "traditional" explanation in terms of dark matter. However, values 
some three to four orders of magnitude larger than that in (^) were needed, 
something which is clearly in extreme disagreement with observations. In 
fact, there seems to be a further mistake, a sign error, in [BJ. A positive 
cosmological constant, as favored by the observations, will tend to accelerate 
the expansion of the universe and, if anything, make matter in the outer 
regions of galaxies less rather than more bound. 

While the effects of a cosmological constant thus are negligible on the 
length scale of galaxies, one might expect observable consequences for galaxy 
clusters. As a first attempt to see such an effect, we will consider the fate of 
circular orbits in a flat, expanding universe with a cosmological constant. To 
do this we start with the equation of motion for a particle in an expanding 
universe with an additional gravitational potential, 
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x + —x = — , 
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where x is the comoving coordinate and a is the scale factor. In terms of 
physical distances, R = a% one finds 



R-R- = -g. 

a 



As an example we consider the de Sitter case with Q\ = 1, i.e. a universe 
totally dominated by the cosmological constant. We then use 



Ht 



a(t) = e 

where H = is constant, to obtain 

v 2 G N m A 

R = ~ 3 K (8) 

for an object in orbit around a central mass m. One might note that the 
same result may be obtained by starting with the static form of the de Sitter 
metric, i.e. 

IG^m A 2 \ 2 dr 2 



ds*=[l- - -r 2 dt l - -P _ - r 2 dn 2 , 

' 1 y _ 2G N m __ A r 2 x 



and using a Newtonian analysis. This form of the metric is related to the 
cosmological form through a coordinate transformation. Equation shows 
that for large enough R, i.e. 



R > 



1/3 



there are no longer bound orbits. Does this have observable consequences? 
Unfortunately it is easy to see that the effect becomes important only for 
orbits that are such that they have periods of the order of the age of the uni- 
verse. Furthermore, the effect of the cosmological constant decreases rapidly 
for smaller orbits. Hence the concept of a rotation curve loses its meaning and 
we had better look elsewhere for a better approach to the possible effect of a 
cosmological constant. One possibility is the way clusters and superclusters 
of galaxies form. 

As a first step, one may consider the infall of matter onto a galaxy cluster 
in the regime where linear perturbation theory is valid. This has in fact 
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been treated by Peebles ||, in the case Qm + = 1 (i.e. zero curvature) 
which is the natural case in view of inflation, and which, incidentally, is now 
also indicated by the recent balloon measurements of the cosmic microwave 
background JF|. Peebles showed that, unfortunately, the dependence of the 
infall peculiar velocity v for a cluster of proper radius R and overdensity 5 
on Q\ is quite weak, being well parametrized by 

v = 0.3H R5Q ^ (9) 

for 0.03 < Qm < 0.3 and 1 < 5 < 3, essentially independent of ^a. This 
formula was generalized to arbitrary Qm + ^a in |§. 

In the non-linear, collapsing phase, we may obtain an estimate of the 
influence of A by adapting the simple constant density, spherical collapse 
model |§ to the situation when the cosmological constant is present. We 
thus look at the situation when an overdense region expands to a maximal 
radius -R max , and then contracts to a viral radius i? v ir- 

To analyse this situation, we use the equation for the energy per unit 
mass (c/. ||, Eq. (20) or our equation (H)) of a mass shell of proper radius 
R(t) containing a fixed mass m: 

„ R 2 G N m AR 2 

E = T"~R 6~~' (10) 

where the three terms correspond to the kinetic energy, Newtonian gravita- 
tional energy, and vacuum energy, respectively. As in the standard analysis 
0, we may employ the virial theorem relating the average value of the kinetic 
energy T to the potential energy V 

dV 

< 2r > = <?• a?>- C 11 ) 

Taking the average over a sphere of constant density of the energy equa- 
tion (at the turn-around radius -R max the kinetic energy is zero) 



E — T v ; r + V vh + V vil — V max + V max , (12) 



utilizing 



and 
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(15) 



we recover in the case A = the well-known result || 

3G N M _ 3G N M 

5-Rmax 10-R v i r 

i.e., i? vir = _R max /2. For a non-vanishing A, the corresponding equation is 

3G N M A 2 3G N M A 
5/? max 10 10.R v i r 5 

Suppose that the mass overdensity contrast compared to the cosmological 
average mass density at the maximal (turn-around) radius is u (= 5.6 in 
the standard case), and assume Qm + ^a = 1. Then we can write M = 
^Rm&x^PM /3, and this inserted into fllTf ) gives, by use of (^) 

1 + k = — + 2/t/i 2 (17) 

where we have introduced \i = R v i r /R max (= 0.5 in the standard case) and 

1 
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This result is written in a somewhat different form than, but agrees with, 
Eq. (26) of ||. If we assume as a first approximation u = 5.6 as in the A = 
case, and = 0.7, we find by solving (|17|) numerically for z = that \i has 
decreased from 0.5 to around 0.39. 

This means that the virialized radius is smaller, which is not unreasonable, 
since for a given mass, only more compact clusters can "survive" the repulsive 
force from a positive cosmological constant. 

We can improve this analysis somewhat by taking into account the fact 
that uj also depends on A. It can be seen that this has the effect of increasing 
uj at intermediate redshifts. This increase causes an increase in fi, meaning 
that we have overestimated the effect of A above. The effect is, however, 
small. One can also estimate what the final density contrast of the cluster 
will be. To do this we have obtained an expression comparing the density at 
maximum expansion with the density of the universe today. This is given by 

~ _ 9vr 2 1 1 

W "T6"/(a) 2 (l-^A) (18) 
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with 



da 



(19) 



(1 - tt A )/a + n A a 2 



The result is a further net relative compression of the clusters due to the 
cosmological constant above the one given by \i. 

Since we have been dealing with an imagined situation where one can 
compare between a "standard" scenario without a cosmological constant, 
and one where A ^ 0, and found only small effects on cluster scales, it seems 
difficult, but maybe not excluded, to draw conclusions about the value of 
in the real universe where observational uncertainties have to be taken into 
account. The effect to search for, is a tendency for virialized clusters to get 
smaller and more overdense for a positive A. 

It seems clear, however, that the effects on galactic scales are extremely 
tiny and negligible for rotation curves. 
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